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Abstract. A triangulated fixed connectivity surface model is investigated by using 
the Monte Carlo simulation technique. In order to have the macroscopic surface tension 
t, the vertices on the one-dimensional boundaries are fixed as the edges (=circles) of 
the tubular surface in the simulations. The size of the tubular surface is chosen such 
that the projected area becomes the regular square of area A. An intrinsic curvature 
energy with a microscopic bending rigidity b is included in the Hamiltonian. We found 
that the model undergoes a first-order transition of surface fluctuations at finite 6, 
where the surface tension r discontinuously changes. The gap of r remains constant at 
£H . the transition point in a certain range of values A/N' at sufficiently large N' , which is 

O | the total number of vertices excluding the fixed vertices on the boundaries. The value 

of t remains almost zero in the wrinkled phase at the transition point while r remains 
negative finite in the smooth phase in that range of A/N'. 
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1. Introduction 

Membranes are known to have a large variety of shapes, which change from one to the 
other depending on some external conditions such as mechanical conditions or some 
chemical environments [U EJ [3] . One interesting topic in membrane physics is to clarify 
how the membrane shape is maintained and to understand the shape of membranes in 
terms of the notion of phase transitions [U [5j [6j [7J . For this reason, the conventional 
surface model of Helfrich and Polyakov 0, El [TUJ has been investigated so far by many 
groups [HI 121 131 [H [15] . Including the conventional model we have currently a variety 
of surface models; extrinsic or intrinsic curvature models p2J [TFJ, [TBI EHJ EDI EH 122] . 
two-dimensional curvature or one-dimensional curvature models [23[ EH [25], fixed 
connectivity models or fluid surface models. 

Phase transitions in the surface models are driven by thermal fluctuations and are 
dependent mainly on the temperature T, because the bending rigidity b[kT], which is 
a parameter included in the Hamiltonian, has the unit of kT. Since b is a microscopic 
parameter, the detailed information of the dependence of b on T is unknown, however, 
b is expected to vary with T. 

The phase transition can also be influenced by the boundary conditions, which 
impose a constraint on the dynamical variable X the position of the surface. The shape 
of membranes depends just on the value of X, and hence the phase structure crucially 
depends on the external conditions that place a constraint on X. A flow field is expected 
to influence the phase structure of the surface model [26J . 

An example of the boundary conditions is to fix a part of the surface. Fixing two 
vertices separated by a distance L on a spherical surface, we have a surface model with 
such nontrivial constraints on X. This model leads us to calculate the macroscopic 
string tension cr of the surface by equating exp(— o~L) ~ Z at sufficiently large L, where 
Z is the partition function of the surface model [271 128]. The expected scaling relation 
of a with respect to N is observed, where N is the total number of vertices of the 
triangulated surface [29l [30j EH [32] . 

It is also possible to fix the tubular surface with one-dimensional boundaries of area 
A. Then, the macroscopic surface tension r is expected to be obtained by exp(— tA) ~ Z 
at sufficiently large A. Therefore, it is interesting to study a surface model that has a 
phase transition and see how the surface tension reflects the phase transition. However, 
r was reported only in [27J. 

In this paper, we study an intrinsic curvature model with two one-dimensional 
boundaries, which are a couple of circles separated by L. Therefore, the surface spans a 
cylindrical tube from one circle to the other. The reason why we use a surface with two 
boundaries is because it seems more convenient to use a surface with two boundaries 
than a surface with one boundary for experimental measurements of the surface tension. 
Let R be the diameter of the boundary circles, then the projected area A is given by RL, 
which is proportional to the area irRL of the tube. In this paper, R is chosen to be R = L, 
then we have A = L 2 . The area nRL of the tube can also be used for the calculation 



Surface tension in an intrinsic curvature model 



3 



of t, however, we use the projected area L 2 for simplicity. We should emphasize that 
it is nontrivial whether the model undergoes a phase transition because the surface is 
fixed by one- dimensional boundaries, which include many vertices, and the size of the 
boundaries increases with the system size, and therefore, the phase structure of the 
model without boundaries [18j or with point boundaries [31] is expected be influenced 
by the one-dimensional boundaries. The purpose of the study is to see whether or 
not the model undergoes phase transitions, and moreover, to see how the transition is 
reflected in the surface tension if the transition occurs in the model. 

2. Model and Mor 
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Figure 1. (a) A rectangular surface of size L\ x L2, and (b) a tubular surface of 
length L\ and diameter L\, where L2 — ~kL\. The projected area of the tube in (b) is 
given by A = L\ . 

Firstly, we show how to construct the triangulated surface on which the model is 
defined. Figure DJa) shows a rectangle of size (Lj, L 2 ). Then we obtain a tubular surface 
in Fig. \Mb) from the rectangle in Fig. []Ja) by bending and sewing up a couple of sides 
of the rectangle. The length L 2 of the rectangle in Fig. [T^a) is chosen to be L 2 = nLi 
so that the diameter R of the cylinder in Fig. H^b) satisfies R — L\. The rectangle in 
Fig. []Ja) and hence the cylindrical surface in Fig. [H(b) can be regularly triangulated, 
and then we have a triangulated open tubular surface with a couple of one-dimensional 
boundaries. We should note that the surface tension r is obtained by a surface with a 
single one-dimensional boundary such as the rectangle in Fig. [T^a) at sufficiently large 
area L\L 2 . However, we should remind ourselves of that it is also possible to define r 
on a tubular surface such as the one in Fig. UJb). 

Figures EJja) and[2^a) show a triangulated surface of size (L 1? L 2 ) = (5, 16) and that 
of (Li, L 2 ) = (30, 94). Because L\ and L 2 are integers, L 2 is given by L 2 = [nLi], where 
[x] denotes the integer obtained by rounding x. The coordination number q is q = 6 
except at the boundary vertices; this can be seen in the lattice in Fig. EJ^a). We should 
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Figure 2. Snapshots of lattices of size (a) N — 80, where (Li, L2) — (5, 16) and (b) 
7V = 2820, where (L x , L 2 ) = (30, 94). 

note that Li = 5 in Fig. [2](a) denotes the total number of vertices along the axis of 
the tube and hence deviates from the length L = 4 by one. The final numerical results 
depend only on the total number of vertices and on the reduced projected area, which 
will be introduced below. 

We multiply Li by a scale factor r in order to vary the projected area A such that 
A = (j-Li) 2 . Table [T] shows the numbers N, N', Lx, L 2 , and the projected areas Ax, 
A2, and A3, and the factors r, which characterize the triangulated tubes used in the 
simulations. N'(=N—2L 2 ) in Tabled] is the total number of internal vertices. The total 
number of vertices N(= LxL 2 ) includes 2L 2 the total number of the boundary vertices. 
Three types of the projected areas Ai(i = 1,2,3) are assumed for the surfaces of size 
iV = 2820 and iV = 5040, and two types of A^i = 1, 2) are assumed for the iV = 9460 and 
iV = 17550 surfaces. The scale factor r iy which is given by 7*j = y/Ai/Li, is as follows: 
n = 17/75 ~ 0.227 and r 2 = 19/75 ~ 0.253 for A = 289 and A = 361 on the N = 17550 
surface, respectively. 

Table 1. The numbers that characterize the triangulated tubes for the simulations; 
the total number of vertices N including the boundary vertices, the total number of 
internal vertices N', the length L\ of the cylinder, the circumference Li of the cylinder, 
and the projected areas A\, A2, and A3, and the scale factors ri = s/Ai/Li. 



N 


N' 


Li 


L 2 


A x 


A 2 


^3 


r x 


r 2 


^3 


2820 


2632 


30 


94 


100 


144 


196 


0.333 


0.4 


0.467 


5040 


4788 


40 


126 


144 


196 


256 


0.3 


0.35 


0.4 


9460 


9116 


55 


172 


196 


256 




0.255 


0.291 




17550 


17082 


75 


234 


289 


361 




0.227 


0.253 





The model is defined by the Hamiltonian, which is a linear combination of the 
Gaussian bond potential Si and the intrinsic curvature energy 5*2 such that 

S = Sx + bS 2 , Sx = J2( X i- X if> S 2 = J2(S i -A l ) 2 , (1) 

(ij) i 
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where <5j is denned by the summation over the vertex angle of the triangles meeting at 
the vertex i. The symbol Aj is defined by 

{2n (i = internal vertex), 
(2) 
7r (i = boundary vertex). 

On the boundaries Aj is fixed to Aj = tt, while it is fixed to Aj = 2ir at the internal 
vertices. 

The partition function Z of the model is defined by 

N' 

n^exp[-S(X)], (3) 

i=i 

where Xj(G R 3 ) is the position of the internal vertex i on the triangulated surface. The 
multiple integrations in Z are performed by leaving the boundary vertices fixed. 
The surface tension r of the model is given by 

2( Sl )-3N' 

2A ' 1 ] 

where A is the projected area of the surface and is given by A— (rLi) 2 denoted as above. 
The expression of r in Eq.(j4]) is obtained by the scale invariant property of the partition 
function [271 125] • 

The canonical Monte Carlo (MC) simulation technique is adopted to simulate the 
multiple integrations WidXi with the Boltzmann weight exp [— S(X)] in the partition 
function. The vertex position X is shifted to a new position X' = X + AX with a 
three-dimensional small random vector AX. The new position X' is accepted with the 
probability Min[l, exp(AS)], where Z\S' = S'(new) — S'(old). 



3. Results 

We show snapshots of surfaces in Figs. [3]Ja) and [3](b) and the surface sections in Figs. 
Efc) and[^d), where the surface size is iV = 17550. The snapshots in Figs. [3](a) and[^c) 
are obtained in the wrinkled phase at 6 = 11, and those in Figs. [3](b) and[3](d) are in the 
smooth phase at 6= 12. Both surfaces are corresponding to the projected area A = 361. 
We find that the surfaces are spanning from one boundary to the other boundary and 
do not shrink to one-dimensional object, and therefore the surface tension r is expected 
to be evaluated in both phases. 

Figures ID^a) - IH(d) show the intrinsic curvature energy S2/N vs. 6, which are 
obtained on the surfaces of size N = 2820, N = 5040, N = 9460, and N = 17550, 
respectively. The corresponding specific heat C$ 2 is defined by 

Cs 2 = ^((S 2 -(S 2 )) 2 ), (5) 

and it is shown in Figs. HJe) -|H(h). We see in Cg a an anomalous peak, which reflects 
a phase transition, where " anomalous peak" denotes that the peak value increases with 
increasing N. We see that the peak value remains low in the case of the largest surface 
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Figure 3. Snapshots of surfaces of size N = 17550 obtained at (a) b— 11 (wrinkled 
phase) and at (b) b— 12 (smooth phase), and (c), (d) the corresponding surface sections. 
The figures are drawn in the same scale. 
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Figure 4. The intrinsic curvature energy S2/N vs. b obtained on the surfaces of size 
(a) N = 2820, (b) N = 5040, (c) N = 9460, and (d) N= 17550, and (e)-(h) the specific 
heat Cs 2 vs. b on the same sized surfaces. The error bars denote the standard errors. 
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of N = 17550 with A = 289. This is because the values of 6, assumed in the simulations, 
are out of the transition region, which is narrow in such a large surface. The transition 
appears to be of first-order, because S2/N discontinuously changes on the surfaces of 
iV = 9460 with A = 196 and iV= 17550 with A = 289 in Figs. 0(c) and 3(d). 

We must comment on the value S2/N, where N includes the boundary vertices, 
on which the definition of 5*2 is slightly different from that on the internal vertices. 
Therefore, S2/N shown in Figs. 0(a) - 0(d) are slightly larger (smaller) than the 
curvature energy per vertex on the boundary (internal) vertices. However, this deviation 
is negligible in the thermodynamic limit, because the ratio L 2 /N', which is the total 
number of vertices on the boundaries over the total number of internal vertices, becomes 
smaller and smaller as the surface size N increases. 

The errors 5s 2 of S2 shown as the error bars in the figures are the so-called the 
standard errors and are defined as follows [33]: Let {S2i}(i = 1, ...,N to t) be a sequence 
of MC data of S 2 obtained at every 1000 MCS, and N tot be the total number of data 
in the sequence. The total number of MCS performed is thus given by 1000xiV tot . The 
series {S 2 i} is split into n b sub-series, and S 2 i denotes the mean value of the 7-th sub 
series. Thus, we have a series of mean values {S2i}{I = 1, n&). Then, 5s 2 is defined by 
the standard deviation of {S 2 i} such that 5$ 2 = \JjyiLi(S2i — S2) 2 / nb, where S2 denotes 
the mean value. In contrast to the standard deviation of the series {S 2 i}, Ss 2 decreases 
if rib is sufficiently large such that {5*2/} are statistically independent. In the analysis 
of data, rib is fixed to «j = 1 x 10 4 in this paper, then the total number of MCS in a 
sub-series is 1 x 10 7 . The error bars on C$ 2 are also the standard errors. 

The two-dimensional bending energy S3 is defined by 

^3 = ^(l-n i -n i ), (6) 

where rij is a unit normal vector of the triangle i, and denotes the summation over 
the nearest neighbor triangles i and j. The surface fluctuations can be reflected in S3. 
Figures 0(a) -0(d) show the two-dimensional bending energy S 3 /N B vs. b obtained 
on the surfaces of size N = 2820, N = 5040, N = 9460, and N = 17550, respectively. 
Nb denotes the total number of internal bonds, where S3 is obtained. The variance 
C*5 3 = ( (S 3 — (Ss)) 2 )/N' is shown in Figs. 0(e) -0(h). We see that the discontinuity in 
S3/NB in Figs. 0(c) and 0(d), and the discontinuity is more accurate than that in the 
curvature energy S2/N in Figs. 0(c) and 0(d). The anomalous peaks seen in C$ 3 are 
also more accurate than those in C$ 2 - The results shown in Figs. 0(a) -0(h) strongly 
indicate that the model undergoes a first-order transition of surface fluctuations. 

The phase transition indicated by a discontinuity of S2/iV and that of S^/Nb is 
expected to be reflected in the surface tension r defined in Eq. (0). Figures 0(a) -0(d) 
show the surface tension r vs. b, which is obtained on the surfaces of size N = 2820, 
iV = 5040, iV = 9460, and iV= 17550, respectively. The surface tension discontinuously 
changes at the transition point b c , where S 2 /iV and S 3 /N B discontinuously change. 
The gap seen in r indicates that the transition is characterized by a difference in the 
macroscopic surface tension. 
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Figure 5. The two-dimensional bending energy S^/Nb vs. b obtained on the surfaces 
of size (a) TV = 2820, (b) iV = 5040, (c) iV = 9460, and (d) N = 17550, and (e)-(h) the 
variance Cs 3 vs. b on the same sized surfaces. 




Figure 6. The surface tension r versus b obtained on the surfaces of size (a) TV = 2820, 
(b) iV = 5040, (c) N = 9460, and (d) N = 17550, and (e)-(h) the Gaussian bond potential 
Si/N vs. b on the same sized surfaces. 
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The surface tension is calculated by using the value of the Gaussian bond potential 
Si, therefore it is expected that Si also discontinuously changes. Figures [6](e) - [6](h) 
show the Gaussian bond potential Si/N' vs. b obtained on the surfaces of size iV = 2820, 
iV = 5040, iV = 9460, and N = 17550, respectively. We see the expected gap of Si/N' in 
Figs. Efg) andEJ^h). We should note that Si/N' — 1.5 is satisfied in the case of a model 
without boundaries such as the spherical surface model because of the scale invariant 
property of the partition function. The reason why Si/N' deviates from Si/N' — 1.5 
shown in Figs. EJ^e) -EJ^h) is because of the constraint on X imposed by the boundaries. 
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Figure 7. Linear plots of r vs. A/N 1 obtained in the smooth phase (Q) and in the 
wrinkled phase (A). 



The surface tension r is expected to scale according to 

r + r = (A/NT, (7) 

where r is a constant that depends on b [271128] . The area A was defined in the previous 
section by using the parameter r such that A— (rLx) 2 , where Li is the distance between 
the boundaries shown in Fig. [IJ Since uLi can be identified with L 2 , then N' is given 
by N = L X L 2 -2L 2 = 7rLi(Li -2) ~ -kL\. Thus we have A/N' = r 2 /ir, and then we 
should vary r in order to plot t+t against A/N'. It is possible to vary r with fixing N', 
however, we vary r by changing N' because we are interested in the size dependence of 
r. 

Since To is unknown, we show in Fig. [7J the surface tension r vs. A/N' in the 
wrinkled phase (O) and in the smooth phase (A). The data are obtained on four 
types of surfaces (N',A) = (17082,289), (N',A) = (9116,196), (N, A) = (4788,144) 
and (N',A) = (2632,100), which correspond to the ratios r = v / 289/75 ~ 0.227, 
r = v / 196/55 ~ 0.255, r = yl44/40 = 0.3, and r = ^/T00/30 ~ 0.33, where r is given 
by r = y/A/Li. The tension r in the smooth phase is obtained at b = 13, while r in 
the smooth phase is at the range 10 < b < 11.5, where the surface is expected to be 
sufficiently wrinkled in each N'. 
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From the results shown in the figure, we find that r remains almost constant both in 
the smooth and in the wrinkled phases in the range 0.016 < A/N' < 0.038. Moreover, r 
remains almost zero in the wrinkled phase while it remains negative finite in the smooth 
phase. Although the phase transition is of first-order, it is not so strong, and then the 
gap of the surface tension is not so clear. For this reason, the results shown in Fig. [7] 
seem not so accurate. Nevertheless, we can conclude from the results that the surface 
tension r discontinuously changes at the transition point, and the value of r in each phase 
remains constant, which is independent of A/N' in the range 0.016 < A/N' < 0.038 at 
sufficiently large surfaces. 



b c (A,N') C 

i — 1 — i — ' — i — 1 — i — 1 — i 




(a) A (b) l-b/b c 

Figure 8. (a) Linear plots of b c (A, N') vs. A obtained on the surface of size N' = 2632, 
iV' = 4788 and jV' = 9116, and (b) log-log plots of C!? ax against l-b/b c . 



Finally, we comment on a scaling property of the specific heat Cs 2 and the strength 
of the phase transition. As we see in Figs. IH(e)-|I](h), the specific heat has the peak 
C!? ax at b c (A, N'), which is the critical bending rigidity and depends on both of the area 
A and the surface size N'. We show in Fig. [5(a) linear plots of b c (A,N') against A, 
which correspond to the data in Figs. 11(e), H](f) and|U(g). The straight lines in Fig. 
Ufa) are drawn by the linear fit of the data, then we have 6 C ~15.5, which is seen to be 
independent of N'. Therefore, 6 C ~15.5 is considered to be the critical bending rigidity 
of the model in the limit of N' — > oo. We should note that the same value of b c can 
also be obtained from Q? ax in Figs. [5(e)-[5(g). By using this value of b c , we plot C^ ax 
against l — b/b c in a log-log scale in Fig. [HJb). The straight lines in Fig. Mjo) were 
drawn by fitting the data to Cl? ax ~ (1 — b/b c )~ u ; the exponent v should be v — 1 if 
the transition is of first-order. We have v = 0.45 ± 0.35 on the N' = 2632 surface and 
^ = 0.76 ± 0.53 on the A^' = 4788 surface. We find that a increases as the surface size N' 
increases, v seems less acurate on the surface of A' = 9116 and is not obtained because 
the total number of data (= 2) is too small for the fitting. Nevertheless, the results 
shown in Fig. [8](b) are considered to be consistent with the first-order transition. 
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4. Summary and conclusions 

In this paper we study a triangulated tubular surface model with a pair of fixed one- 
dimensional boundaries by using the canonical Monte Carlo simulation technique. The 
model is defined by the Hamiltonian that is a linear combination of the Gaussian bond 
potential Si and the intrinsic curvature energy S2. The model is known to have first- 
order transitions on the surfaces without boundaries and on those with point boundaries 
|18L [3Tj . However, it is non-trivial whether the model undergoes a phase transition with 
the one-dimensional boundaries. The surface shape of the model in this paper is almost 
fixed by the one-dimensional boundaries, and hence, the transition of surface fluctuations 
is expected to be influenced by such boundary conditions because the transition is always 
accompanied by the collapsing transition, which is a transition of shape transformations. 
Therefore, special attentions are paid on whether the model undergoes a phase transition 
and how the phase structure is reflected in the macroscopic surface tension r. In order 
to compute r, we assume several values of the projected area A of the tubular surface of 
a given size N'. The surface tension is expected to scale according to r + r = (A/N') u 
with a constant tq in the limit of N' —>■ 00 and A — > 00 while A/N' is fixed [271 28J. 
However, we are interested in the dependence of r on A/N' rather than the scaling 
behavior because no information of r is obtained. The assumed values of A/N' are in 
the range 0.016 < A/N' < 0.038. We use the surfaces of size up to N' = 17082, which 
excludes the total number vertices on the fixed boundaries. 

Firstly, we find that the model undergoes a first-order transition of surface 
fluctuations, which separates the smooth phase from the wrinkled phase. Both of 
the potential Si and the intrinsic curvature energy S 2 discontinuously change at the 
transition point b c on relatively large sized surfaces in the range 0.016 < A/N' < 0.038. 
The two-dimensional bending energy S3, which is not included in the Hamiltonian, has 
also a gap at the transition point b c , and this indicates that the smoothness of the surface 
discontinuously changes at b c . Secondly, the surface tension r is also found to change 
discontinuously at the transition point b c . The value of r is constant in the smooth 
phase and it is almost zero in the wrinkled phase at 6(7^ b c ) close to the transition 
point in the range 0.016 < A/N' < 0.038 at sufficiently large N', and hence the gap of 
t remains finite at the transition point. Finally, we find by finite-size scaling anaysis 
that the specific heat C!? ax scales according to an expected relation. The results are 
consistent with the first-order transition. 

It is interesting to study whether the transition is seen and reflected in the surface 
tension in the same model on fluid surfaces. It is also interesting to study connections 
between the surface tension and the phase transitions in the other type of surface models. 
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